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One-dimensional quasiperiodic systems, such as the Harper model and the Fibonacci quasicrystal, 
have long been the focus of extensive theoretical and experimental research. Recently, the Harper 
model was found to be topologically nontrivial. Here, we derive a general model that embodies a 
continuous deformation between these seemingly unrelated models. We show that this deformation 
does not close any bulk gaps, and thus prove that these models are in fact topologically equivalent. 
Remarkably, they are equivalent regardless of whether the quasiperiodicity appears as an on-site or 
hopping modulation. This proves that these different models share the same boundary phenomena 
and explains past measurements. We generalize this equivalence to any Fibonacci-like quasicrystal, 
i.e. a cut and project in any irrational angle. 
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Recent experimental developments in photonic crys- 
tals [H [5] and ultracold atoms [3HS] have made the study 
of the dynamics of particles in one-dimensional (ID) 
quasiperiodic systems experimentally accessible. These 
fascinating systems have long been the focus of extensive 
research. They have been studied mainly in the context 
of their transport and localization properties, showing a 
variety of interesting transitions between metallic, local- 
ized, and critical phases [MSj . With their recently found 
nontrivial connection to topological phases of matter [TO] 
there is a growing interest in their boundary phenom- 
ena QHCEU. 

The behavior of particles in such systems is described 
by ID tight-binding models with quasiperiodic modula- 
tions. There is an abundance of quasiperiodic modula- 
tions, among which the canonical types are the Harper 
model (also known as the Aubry- Andre model) [SJ [T3] 
and the Fibonacci quasicrystal (QC) |14) . The quasiperi- 
odicity of the Harper model enters in the form of a co- 
sine modulation incommensurate with lattice spacing, 
whereas the Fibonacci QC has two discrete values that 
appear interchangeably according to the Fibonacci se- 
quence. Moreover, the quasiperiodicity may appear in 
on-site terms (diagonal), in hopping terms (off diago- 
nal), or in both (generalized). Each of these models de- 
scribes different physical phenomena. Indeed, the Harper 
and the Fibonacci modulations have different localization 
phase diagrams, depending also on their appearance as 
diagonal or off-diagonal terms (see e.g., Refs. [7H51 \TE\). 
Notably, several attempts were made to gather these 
models under some general framework [9l I16[ 117) . but 
with only partial success. 

These ID quasiperiodic models play a nontrivial role 
in the rapidly growing field of topological phases of mat- 
ter [T0H12] . This new paradigm classifies gapped systems 
such as band insulators and superconductors [HI IT9] . 
Each gap in these systems is attributed an index that 
characterizes topological properties of the wave functions 
in the bands below this gap. By definition, two gapped 
systems belong to the same topological phase if they can 



be deformed continuously from one into the other with- 
out closing the energy gap. Conversely, while deforming 
a system with a given topological index to a system with 
another index the bands invert and the bulk gap closes, 
i.e. a quantum phase transition occurs. 

The spectra of the aforementioned ID quasiperiodic 
models are gapped and, hence, appropriate for topolog- 
ical classification. Apparently, in the absence of any 
symmetry all ID systems are topologically trivial [20) : 
namely, all their topological indices are zero, and there- 
fore the hopping terms in such systems can be contin- 
uously turned-off (the atomic limit). Conversely, under 
the same conditions 2D systems have nontrivial topologi- 
cal phases, which are characterized by an integer index - 
the Chern number. If a Hamiltonian of a ID system de- 
pends on a periodic parameter, then this parameter can 
be considered as an additional dimension. Taking into 
account all the possible values of this parameter, the sys- 
tem becomes effectively 2D, and may have a nontrivial 
Chern number [21] , This is seen during the evaluation 
of the Chern number, which requires integration of the 
Berry curvature over this parameter. It was shown in 
Ref. [TU] that for QCs, a translation of the quasiperiodic 
modulation may be considered as such a parameter. Re- 
markably, for QCs the Berry curvature is invariant to 
this translation, making the integration over it redun- 
dant. Therefore, it was concluded that ID quasiperiodic 
systems can be associated with Chern numbers, making 
the ID QCs topologically classified. 

Here we use this approach to prove that all the afore- 
mentioned ID models are topologically equivalent when- 
ever they have the same modulation frequency. Hence, 
there is no quantum phase transition when deforming be- 
tween these different models. To prove this equivalence, 
we extend each one of the ID models to an "ancestor" 2D 
model. We find that all the resulting models are variants 
of the 2D integer quantum Hall effect on a lattice. These 
2D models are topologically equivalent and nontrivial. 
Therefore, their corresponding ID descendants are also 
topologically equivalent and nontrivial. Remarkably, the 
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equivalence holds between any Fibonacci-like quasicrys- 
tal and a Harper model with a corresponding modulation 
frequency. 

One-dimensional tight-binding Hamiltonians with 
nearest-neighbor hopping and an on-site potential can 
be written in the general form 

H = E [ (* + A ° d ^° d ) c » c "+i + h - c - + A X d 4c„] , (i) 

n 

where c„ is the single-particle annihilation operator at 
site n, t is some real hopping amplitude, V od is some hop- 
ping modulation (off-diagonal term) and V d is an on-site 
potential (diagonal term). The real and positive param- 
eters A od and A d control the strength of the off-diagonal 
and diagonal modulations, respectively. The quasiperi- 
odicity of the different models is encoded in the poten- 
tial modulations, V od and V d . Using this general form, 
we first show the topological equivalence between all the 
Harper models, since their relation to the quantum Hall 
effect is evident. Then we turn to encompass also the 
Fibonacci QCs. 

Let us begin with the Harper models, which are gov- 
erned by the modulation V^(k) = cos(2irbn + k). This 
modulation is parametrized by the frequency b and phase 
k. Whenever b is irrational, the modulation is incommen- 
surate with the lattice and describes a QC. In this case, k 
resembles a translation of the quasiperiodic modulation. 
Since k does not affect bulk properties, it was usually 
ignored in previous analyses. However, as we shall soon 
observe, it plays a crucial role in unraveling the topolog- 
ical behavior of quasiperiodic models [TU] . 

The diagonal Harper model 1 1 3 j is defined by setting 
A od = 0, A d ^ 0, and V d = V^{k) in Eq. 0. The cor- 
responding Hamiltonian describes uniform hopping and 
a modulated on-site potential. For any given k, this 
Hamiltonian can be viewed as the A:th Fourier compo- 
nent of some ancestor 2D Hamiltonian. From this view- 
point, k is a second degree of freedom, hence we define 
the operator c n ^ that obeys the commutation relation 
{ c n,k,c n i k'} = $n,n'3k,k'- We can now define a 2D Hamil- 
tonian % = /p 71 (dk/2ir)H(k), where in H(k) we replace 
the operators c n with c n ^. Note that, in the following, H 
(H) denotes ID (2D) Hamiltonians. Defining the Fourier 
transform c„.fc = J2 m e ~~ lknic n,m, we obtain the 2D an- 
cestor Hamiltonian of the diagonal Harper model 

K=Y1 [ tC lrnCn+l, m + ye^Wl + h ' C ] • ( 2 ) 

This Hamiltonian describes electrons hopping on a 2D 
rectangular lattice in the presence of a uniform per- 
pendicular magnetic field with b flux quantum per unit 
cell [H [21 [23], as illustrated in Fig.[jja). Note that, in 
the magnetic field appears in Landau gauge. 
In the absence of a magnetic field, the Hamiltonian 
commutes with the group of translations. Since the 




FIG. 1. Graphical presentations of 2D Hamiltonians. The 
electrons hop on a rectangular lattice in the presence of a 
perpendicular magnetic field with b flux quantum per unit 
cell. A hopping amplitude from one vertex to another is de- 
noted by an arrow, (a) The Hamiltonian the "ancestor" 
of the diagonal Harper model. The hopping is to nearest 
neighbors, and each rectangular plaquette (marked in cyan) 
is pierced by b flux quantum, (b) The 2D ancestor Hamilto- 
nian of the off-diagonal Harper model, "H^ d . The hopping is to 
nearest neighbors in the n direction, and also to next-nearest 
neighbors. Here, there are two types of plaquettes, a triangle 
(dark green) and a parallelogram (cyan), both pierced by b 
flux quantum. The magnetic translation group is the same 
in both Hamiltonians, implying that they are topologically 
equivalent. 

magnetic field breaks this symmetry, the notion of the 
magnetic translation group was introduced J2jT. The 
magnetic translation group is generated by the oper- 
ators Tfn and T A , where T^Cn^T^ 1 = c n , m+ i and 
ThCn im T7 l = e~ t2nbm c n+ i jm . These operators commute 
with the Hamiltonian but not with each other. However, 
for a rational flux, 6 = p/q, the operator T q n = {Tn) q 
commutes with T^. Therefore it is possible to diago- 
nalize simultaneously , T q n, and T^. The spectrum 
in this case is composed of q bands [23 . In a seminal 
paper by Thouless et al. [25J, it was shown that each 
gap in the spectrum of this model is associated with a 
quantized and nontrivial Chern number (Hall conduc- 
tance). Later on, it was shown by Dana et al. [26] that 
the nontriviality of the Chern numbers stems from the 
symmetry of this model with respect to the magnetic 
translation group. Choosing consistent boundary condi- 
tions [25], the Chern number v r that is associated with 
a gap number r = 1, . . . , (q — 1) abides the Diophantine 
equation r = v r q + t r p, where v r and t r are integers, and 
< \v r \ < q/2 [25]. 

An irrational b can be approached by taking an ap- 
propriate rational limit with p, q — > oo. In this limit the 
spectrum becomes fractal [23]. Nevertheless, even for an 
arbitrarily large q, the system abides the aforementioned 
Diophantine equation, and the gaps remain associated 
with nontrivial Chern numbers. Formally, the evaluation 
of the Chern numbers requires integration of the Berry 
curvature over k [25J. Hence, the quantized Chern num- 
bers characterize only the 2D ancestor model, rather than 
its ID descendant model. However, as shown in Ref. [10j . 
for a QC, i.e., for an irrational 6, the Berry curvature is 
independent of k and there is no need for such an integra- 
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tion. Therefore, the ID models can be associated with 
the same quantized topological indices. 

This simple model demonstrates that the method to 
extract the topological indices of a ID QC is to extend 
it to 2D using the above procedure and find the mag- 
netic translation group of the ancestor Hamiltonian. This 
yields a Diophantine equation and, thus, the Chern num- 
ber of each gap. 

Having performed this method for the diagonal Harper 
model, we turn, now, to the off-diagonal Harper model. 
In this model the hopping is modulated and the on-site 
potential vanishes. It is defined by setting A od 7^ 0, A d = 
0, and V° d = V^(k) in Eq. 0. Constructing its 2D 
ancestor model, we obtain 



r \ od / 

njH _ 1 , t 1 (J2nbnJ 

od / j Y' c n,m c n+l,m 2 V c n,m c n+l,m+l 



— ilixbn f \ ■ r 



(3) 



This Hamiltonian describes electrons hopping on a rect- 
angular lattice with nearest-neighbors hopping only in 
the n direction, and also next-nearest-neighbors hopping. 
Here too a perpendicular magnetic field is present with 
b flux quanta per unit cell and per plaquette, as illus- 
trated in Fig. [ljb). The corresponding magnetic trans- 
lation group is the same as in the diagonal case. Hence, 
this model abides the same Diophantine equation, which 
characterizes its gaps. Therefore, for a given b, the 2D 
ancestor Hamiltonians of the diagonal and off-diagonal 
Harper models have in fact the same number of gaps and 
the same distribution of Chern numbers, making them 
topologically equivalent. 

The diagonal and off-diagonal Harper models are in- 
corporated in the generalized Harper model 8j, where 
in Eq. ([!) we take A od ^ 0, A d ^ 0, V d = V*(k), and 
V° d = V^-{k + irb). Now both the hopping terms and the 
on-site potential are cosine modulated. Its corresponding 
2D model has both nearest-neighbor and next-nearest- 
neighbor hopping. The magnetic flux per unit cell is 
still b, and the magnetic translation group remains the 
same as well. Hence, the Diophantine equation is also 
the same, independent of relative modulation strengths 
A d and A od . We can therefore conclude that, for a given 
b, as the ratio A od /A d is changed all the energy gaps re- 
main open and thus keep their Chern numbers fixed |27j . 
This means that the generalized Harper model provides 
a way to continuously deform the diagonal Harper model 
into the off-diagonal one, and vice versa, without experi- 
encing a quantum phase transition. 

After showing the topological equivalence between the 
Harper models, we now address the Fibonacci QC. This 
QC is governed by the modulation V£ = 2([(n + 2)/rJ — 
L(n + l)/rj) - 1 = ±1, where r = (1 + V5)/2 is the 
golden ratio and |_a;J is the floor function. Similar to the 
Harper modulation, the V F modulation can be employed 




FIG. 2. (a) The smooth continuation between the Harper 
and the Fibonacci modulations [cf. Eq. Q] for b = 
2/(1 + \/5), k — 0, and several values of the smoothing pa- 
rameter ft. For ft — > 0, V s resembles the Harper modulation 
F H , whereas for ft — > oo, it approaches the Fibonacci modu- 
lation V F . Varying ft has no influence on the distribution of 
the Chern numbers, meaning that the Harper model and the 
Fibonacci QC are topologically equivalent, (b) Illustration of 
the space of quasiperiodic systems spanned by the smoothing 
parameter ft and the diagonal-off-diagonal ratio A od /A d , when 
substituted in the general Hamiltonian H s [cf. Eq. |5|]. All 
the systems in this space are topologically equivalent, with 
the same number of gaps and the same Chern numbers. 

as a diagonal Fibonacci QC, with A od = and V d = V* . 
Alternatively, it can be employed as an off- diagonal Fi- 
bonacci QC, with A d = and V° d = V* . Due to the 
discontinuous nature of these QCs, they have no appar- 
ent ancestor 2D models. This seemingly prevents the 
extraction of their topological indices. 

We can overcome this barrier by constructing a mod- 
ulation that continuously deforms the Fibonacci into 
a Harper modulation. Consider the function f(x) — 
2([x + a\ — [x\) — 1, where < a < 1. The function 
h(x) — cos(27ra; + air) — cos(a7r) has the same sign as 
f(x) for any x. Therefore g(x; f3) — tanh[/3/i(x)]/ tanh[/3] 
is a continuation between the smooth h{x) = g(x; ft — > 0) 
and the steplike f(x) = g{x;(3 — > oo). 

Accordingly, we define the smooth modulation 



tanh{/3[cos(27r&n + k) — cos7r6]} 
tanh j3 



(4) 



It can be seen that, in the limit of small /3, this smooth 
modulation becomes the Harper modulation, up to a con- 
stant shift, V£(k;fi -> 0) = V^(k) - cosTrfo. In the 
opposite limit, it approaches the Fibonacci modulation, 
V„(k = 3-7r6;/3 —> oo) = V r f for b — 1/r, as depicted in 
Fig. [2ja) . Similar to the generalized Harper model, we 
now define a generalized smooth model 

H s (k; p) = Y,{ [t + A od F„ S (k + 4tt6; P)]clc n+l + h.c. 
n 

+ \ d V*{k + Sub; p)clc n }, (5) 

where b — 1/t. This model is a continuous deformation 
between the generalized Harper model and a generalized 
Fibonacci QC, with (3 being the deformation parameter. 

Now we are able to construct a corresponding 2D 
model that includes also the Fibonacci QC. Recall 
that, in the diagonal Harper model, the modulation of 
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cos(27r£>n + k) resulted in hopping one site in the m di- 
rection, accompanied with a phase factor of e vbn [see 
Eq. pb]. One can think, however, of a more general mod- 
ulation F[cos(2irbn + k)], where F is some analytic func- 
tion. In such cases, F can be expanded as a Taylor series 
in powers of cos(27rfon + fc). In turn, the series can be 
rewritten as a series in powers of e ± l ( 27rfm + fe ) p The Zth 
component of this series will appear in the 2D ancestor 
model as a term of hopping I sites in the m direction. The 
amplitude of this hopping term would be Fie %2 ' K ° n , where 
Fi is the Fourier transform of F[cos(fc)]. Similarly, an off- 
diagonal modulation will result in terms that incorporate 
hopping to the nearest-neighbor in the n direction with 
longer-range hopping in the m direction. Therefore, the 
2D ancestor Hamiltonian of H s (k;(3) is 



+ l,m 



(6) 



OG 

E^)( A ° d 



iKb(2n+±)l f 



1=0 



+ A od e - fe 6(2n + 4)i c t >jnCB+1(m _ J 

+ X d e^ 2n+ ^ l ci >mCn , m+l ) + h.c. 



where vf(/3) = (dk/2n)e ilk V^ =0 (k; f3). It can be 



shown that 



l e 2Z 
i e 2Z + l 



e tanh P ^J=0 



S tanh I) 



, i+i 



(7) 



where Zj(J$) — yT-j- — ia) 2 — (u>j — ia) with 
u!j = (j + 1/2)tt/(3 and a = cos(7r6). The physical mean- 
ing of vf is better understood by looking at the limits of 
large f3 and small f3. In the Harper limit of /? <C 1 [2S], 



wf()8<l) 



/ G 2Z 
/ G 2Z 



-a 



(8) 



where u; = 1 + (£;,o + <^,±i) ( 7r2 /8 — l) ■ We can see that, 
in this limit, the hopping in the m direction decays ex- 
ponentially with the distance I. In the extreme limit of 
/3 = 0, only the terms with I — 0, ±1 survive and T-L be- 
comes the 2D ancestor of the generalized Harper model 
(up to a constant shift of the energy). In the opposite 
limit of f3 ^> 1, i.e., the Fibonacci limit [2"8] , 



«f(/3 »!)«(-!) 



frin(j7r6)+<Ji 1 o(l + 2|&J) . (9) 



Now the hopping in the m direction is no longer local, 
but decays as 1/Z. 

Regardless of the exact value of (3, W s describes elec- 
trons hopping on a rectangular lattice. Moreover, since 



the amplitudes of the hopping in the m directions are 
accompanied by the phases e ±l27rbn ', in this model also 
a magnetic field is present with b flux quanta per unit 
cell. Remarkably, despite the varying hopping behavior, 
the magnetic translation group remains the same for all 
values of j3. Therefore, the Diophantine equation is also 
the same for all /?. Consequently, while varying (3 from 
zero to infinity the gap structure and its corresponding 
nontrivial Chern numbers are unchanged. Since j3 turns 
the Harper model into the Fibonacci QC, it implies that 
they are topologically equivalent. Note that this is true 
also for any value of A od /A d . 

The real and positive parameters, j3 and A od /A d , span 
the space of models that contain the diagonal and off- 
diagonal, Harper and Fibonacci models, as illustrated in 
Fig. [2jb). We can therefore conclude that all the ID 
models H s in this space are topologically equivalent and 
nontrivial. Moreover, the same holds true for any irra- 
tional b. Here, taking j3 to infinity results in a Fibonacci- 
like QC with t = 1/6 as the modulation frequency of V F . 
Note that any Fibonacci-like QC can be obtained via the 
cut-and-project method 29J, by taking r = 1 + 1/ cot 6 
with 9 the angle of the projection line. This means that 
any Fibonacci-like QC with a given r is topologically 
equivalent to a Harper model with b = 1/r, and both 
are topologically nontrivial. For a rational b, i.e., for pe- 
riodic systems, the 2D ancestor Hamiltonians are equiv- 
alent, but the implications to the ID models are more 
subtle 110 • Nevertheless, no bulk gaps will close when 
deforming between them. 

The physical manifestation of the topological nontrivi- 
ality of a QC is easily seen in the emergence of boundary 
states that traverse the energy gaps as a function of the 
translation parameter k. The Chern number of each gap 
is equal to the number of boundary states that traverse 
the gap. The topological equivalence of the aforemen- 
tioned ID models implies that the number of traversing 
boundary states in a given gap is constant when j3 and 
A od /A d are changed. This should be contrasted to the 
localization properties of the bulk states, which vary con- 
siderably during such a change [8]. 

The topological nontriviality of the diagonal and off- 
diagonal Harper models has been recently demonstrated 
experimentally via boundary phenomena 110] . It would 
be intriguing to test this for the Fibonacci QC as well, 
where we expect gap-traversing boundary modes to ap- 
pear. Our prediction is supported by the fact that the 
existence of subgap boundary states in the Fibonacci QC 
was noticed [3D] , measured [3TJ [35] , and analyzed to be 
quantitatively similar to those of the Harper model [33]. 

To summarize, we developed a ID model that ranges 
smoothly from the Harper model to a Fibonacci or 
Fibonacci-like QC and from diagonal to off-diagonal 
modulations, as a function of control parameters. Us- 
ing the fact that dimensional extension of QCs from one 
to two dimensions reveals their topological character, we 
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extended this model to 2D. We found that in 2D the hop- 
ping behavior changes significantly with the control pa- 
rameters. Nevertheless, the magnetic translation group is 
unaffected. This implies that the same nontrivial Chern 
numbers remain for all values of the control parameters. 
Therefore we conclude that all these ID models are topo- 
logically equivalent and nontrivial. It would be interest- 
ing to follow the same process for other QCs, which may 
exhibit novel types of topological phases. For example, 
2D and 3D QCs which are obtained via cut-and-project 
methods [29, 34], may have topological characteristics of 
4D and 6D systems, respectively. 
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We thank the Minerva Foundation, the US-Israel Bina- 
tional Science Foundation, and the IMOS Israel-Korea 
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SUPPLEMENTARY MATERIAL 



In the main text we defined a quasiperiodic modu- 
lation, V^(k;(3) [cf. Eq. (4) of the main text], which 
smoothly deforms the Harper modulation into the Fi- 
bonacci modulation as j3 is changed from zero to infinity. 
Placing it in the generalized smooth model [cf. Eq. (5) of 
the main text] and constructing an ancestor 2D Hamil- 
tonian requires the evaluation of the following integral 



Ti = 



2tt 



dk 
2t7 



e ilk V^(k + 3Trb; 



(10) 



Recall that V^(k + 3tt6;/3) = tanh{/3[cos(7„ + k) - 
a]}/ tanh /3, where 7„ = 2-7r6n + 37r6 and a = cos irb. 
Shifting the integration variable k — > k — j n and tak- 
ing into account the 2ir periodicity of the integrand, we 
obtain 



Ti = 



tanh/3 



2tt 



dk 
— ( 
2tt 



, i I k 



tanh{/3[cos(fc) — a}} 



(11) 



Since the integrand is real and symmetric, we know that 



>-l 



Therefore it is sufficient to evaluate vf for 



I > 0. Defining z = e 

v s = 1 / — 
1 tanh/3 J unit 2m 



tanh 



2a) 



(12) 



where the integration is performed over the unit circle in 
the complex plane. 

For 1 = 0, there is a pole at z = 0. Below we show 
that this pole has zero contribution. For any I > 0, the 
integrand has simple poles when 



z + z 1 = 2i(wj — ia), 



(13) 



with Wj = (j + l/2)n//3 and j is an integer. Specifically, 
these poles are 



(Wj - ia) ± yj (wj - iaf 



+ 1 



(14) 



The poles that contribute to the integral are only those 
that are within the unit circle. Below, we prove that 

< 1 for j > 0. 

the relevant 



„+l 



: , I < 1 for j < —1, and that 



Therefore, using the fact that z_ ■ = (zj_ 1 
poles are z~ and \zj) for j > 0. 

Denoting zj = —izj, where Zj = (wj — ia) 2 + 1 — 
(wj — ia) , the residue of the integrand at z~ is given by 



Res <! z l 1 tanh 



(z + z^ 1 - 2a) 



Z 3 
i+1 



/3 1 + z. 



(15) 



We therefore conclude that 
a 2i l ~ l 



-7 



/? tanh /3 



E 

3=0 



i+1 



(-1)' 



l + z 2 + l + {z*f 



. (16) 



which is equivalent to Eq. (7) in the main text. 

Having the exact serial expression for vf(/3), we turn 
to approximate it in the limits of small and large j3, i.e. 
the Harper and Fibonacci limits, respectively. This way, 
we obtain simpler expressions, which provide better un- 
derstanding of the physical effect of (3 on the 2D ancestor 
Hamiltonian. 

Small (3 limit (Harper limit) 

In the limit of j3 -C 1, by definition, Wj 3> 1 for every 
j > 0. Therefore, 

^<1) = ^(1 + ^+0(0. W 



and 



i+1 



1 + z 2 (2w j y+ 1 



!+»(/ + !)- 



(18) 



Plugging this into Eq. ( 16 ), we consider the cases of even 
and odd values of I. 

For odd values of I, we obtain 



2-1 oc 



4 2Z+1 (/3« i)^(l)" E^y+r («) 

Notice that E^=o(i+ 1 / 2 )"^ 1 = (2 Z+1 — 1)C(^ + 1) , where 
C{%) is the Riemann zeta function. Using the facts that 
C(2) = tt 2 /6, and that for I > 3, + 1) ~ 1 and 2 l+1 > 
1, we can write 



1*1 = 1 
|J|=3,5, 



«l62Z+lG8« 1) ? 

For even values of I, we obtain 



1/2 

^(ifi) 1 " -1 ■ ( 2 °) 



f 6M 0<l) «-*(/ + !) A (i) l tj7^ W ( 21 ) 



^ 2 w ^(i + i/2) ; + 2 

Similarly, E^oO' + V 2 )^ 2 = (2 l + 2 - l)((l + 2), result- 
ing in 



«fe2z(/3 « 1) 



I = 



|/|=2,4,... -a^O + 1) f 



l<l . (22) 



The results of Eq. ( 20 1 and Eq. ( 22 ) are summarized in 



Eq. (8) of the main text. 

Large (3 limit (Fibonacci limit) 

In the opposite limit of /3 ^> 1, the poles Wj becomes 



7 



dense. Therefore we can approximate the sum YljLo 
with the integral ((3/ it) J„/2p v ?j- Additionally, using the 

1/2 

relation 1 + Zj 2 = 2Zj [l + (wj — ia) 2 ] , we find that 
z /w = -2z - 2 (1 + z 2 )- 1 . Hence, for I > 0, 



oo l + l 



i=0 -+v 



2tt 



l-i ?j 



tt/2/3 



2^ V " ! 



z-'(w 3 = tt/2/3). 



(23) 



In the limit of /? — > oo, we observe that itfj <C 1 for any 
finite I. Since |a| < 1, we obtain 



+ Ok 2 )- 



3=0 



Rewriting y/l + a 2 + ia = e 

00 z 

y — 

For Z = 0, 

oo 

1 + z-- 



% arcsm a 



, we have 



' p— fr-oo P it arcsin a 

z 2 2ttZ 



(24) 



(25) 



3=0 



2tt 



/2/3 



= — [asinh(w)|„,_i. 00 — asinh(7r/2/3 — ia)] 

27T 



2tt 



[asinh(?x;)| l 



i arcsin a 



(26) 



Therefore, according to Eq. (16), we obtain 



Recalling that \a\ < 1, we can bound |i> e | < e 2 , meaning 
that v e vanishes when e —> 0. In this limit the only pole 
of the integrand that resides within C e is z — 0. We can 
therefore conclude that the contribution of this pole to 



the integral in Eq. ( 12 ) is zero. 



Poles within the unit circle 

In this last part we prove that the poles 2j that reside 
within the unit circle are those with j < — 1 for z^ and 
with j > for zj . 

We do so by proving that \zf\ = 1 iff Wj = 0. Note 
that Wj =/= for any integer j, but changes its sign from 
negative to positive only upon an increment from j = — 1 
to j = 0. From Eq. ( 13 1 we can deduce that — > oo 



0, while Za 



and Za 



— > oo. 



and Zj > x , x 

Therefore, this proves the preposition. 

For brevity, we denote u) — wj — ia, and rewrite 
Eq. (13) as z(z — 2iu) = — 1. Multiplying each part of 



this equation by its complex conjugate gives 
\z\ 2 [\z\ 2 + 4M 2 + 4hn(z*w)] = 1. 



(30) 



Requiring \z\ = 1, we obtain \lo\ 2 = -Im(z*w). Since 



z to 



\w\ z* /<*)*, we get lm(z/u>) = 1 



Recalling from Eq. (14) that z — i(cj ± vT+ uj2 ), we 
end up with Re(\/1 + w -2 ) = 0. This is true only if w is 
purely real and u~ 2 < — 1. By definition this is the case 
only when \a\ < 1 (which is consistent) and = 0. 



vf f/3 — > oo) = — sin(Z arccos a) — <5; o, (27) 

where we used the facts that arcsin a = tt/2 — arccosa, 
and that sm(lx)/l = x for I = 0. Recalling that a — 
cosnb, we obtain Eq. (9) of the main text. 



The pole at z = 

We have mentioned above that the integrand in Eq. ( 12 ) 



has a pole at z — for I = 0, and claimed that it has 
null contribution to the integral. In order to prove this, 
we consider the integral 



1 



tanh f3 J c 2ni 



dz 



-z* tanh 



(z + z* - 2a) 



(28) 



where C e is a circle in the complex plane defined by \z\ 
e. By denoting z — ez, we obtain 



dz 

tanh/3 J unit 2iri 
1 '*« dk 



z 1 tanh 



(ez + ez 1 - 2a) 



tanh (3 J 2tt 



tanh{/3[ecos(fc) — a]}. 



(29) 



